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Spherical Relativistic Radiation Flows with Variable Eddington Factor 
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Abstract 

We solve spherically symmetric radiation flows under full special relativity with the help of a variable 
Eddington factor /(r, /3), where r is the optical depth and /? is the flow velocity normalized by the speed 
of light. Relativistic radiation hydrodynamics under the moment formalism has several complex prob- 
lems, such as a closure relation. Conventional moment equations closed with the traditional Eddington 
approximation in the comoving frame have singularity, beyond which the flow cannot be accelerated. In 
order to avoid such a pathological behavior inherent in the relativistic moment formalism, we propose a 
variable Eddington factor, which depends on the flow velocity as well as the optical depth, for the case of 
the sperically symmertic one-dimensional flow. We then calculate the relaticistic spherical flow with such 
variable Eddington factors to investigate the case that gas is accelerated by radiative force. As a result, it 
is shown that the gas speed reaches around the speed of light by radiation pressure. 
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1. Introduction 

Relativistic outflows from a luminous central object are 
observed in various active phenomena; e.g., relativistic 
jets and winds from microquasars (Mirabel, Rodriguez 
1999; Fender et al. 2004), such as SS 433, GRS 1915+105, 
GRO J1655— 40, jets in active galactic nuclei, such as 
3C 273 and gamma-ray bursts (Meszaros 2002). Several 
mechanisms were proposed to explain these relativistic 
outflows, including hydrodynamical, radiative, and mag- 
netic drives. When the luminosity highly exceeds the 
Eddington one, the relativistic outflow seems to be driven 
by radiation pressure of the central object. 

So far relativistic outflows or winds driven by radia- 
tion pressure in the spherically symmetric case have been 
studied by several researchers (Castor 1972; Ruggles, Bath 
1979; Mihalas 1980; Quinn, Paczyhski 1985; TuroUa et al. 
1986; Paczynski 1990; King, Pounds 2003). As for numer- 
ical calculation, the radiation transfer has been sovled in 
two or three dimensions by a Newtonian treatment, but 
has not been resolved sufficiently for a highly relativistic 
case yet. 

On the other hand, under the traditional moment for- 
malism, the relativistic outffows driven by radiation pres- 
sure have pathological behavior (e.g., TuroUa, Nobili 1988; 
NobiH et al. 1991; Turolla et al. 1995; DuUemond 1999; 
Fukue 2005). That is to say, moment equations for rel- 
ativistic radiation transfer can have unphysical critical 
points. For example, in one-dimensional relativistic ra- 
diation ffow using the Eddington approximation in the 
comoving frame, where the moment equations are trun- 
cated at the second order, the singularity appears when 
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the ffow velocity becomes c/^/3. This is understood as fol- 
lows (TuroUa, Nobih 1988; Nobih et al. 1991; DuUemond 
1999; Fukue 2006). The radiative diffusion may become 
anisotropic even in the comoving frame of the gas as a 
result of what the velocity gradient becomes very large in 
the direction of the ffow when the gaseous ffow is radia- 
tively accelerated up to the relativistic regime. Hence, in 
a flow that is accelerated from subrelativistic to relativis- 
tic regimes, the Eddington factor should be different from 
1/3 even in the optically thick diffusion limit. 

As already stressed in the literature (e.g., Nobili et al. 
1991), the location of the critical point in the moment 
equations depends on the choice of a closure relation, and 
with a suitable choice of the closure relation, the critical 
point may disappear. For example, Nobili et al. (1991) 
adopted a variable Eddington factor, which depends on 
the optical depth. However, for the present transfer ffow, 
the critical condition, where the denominator of moment 
equations vanishes, contains the flow velocity (e.g., Nobili 
et al. 1991; Fukue 2006). Hence, as a natural extension, in 
the present study we have proposed a variable Eddington 
factor which depends on the "flow velocity" as well as 
the optical depth. By adopting such a velocity-dependent 
Eddington factor, we intend to send the critical point away 
toward the edge of the speed of light. 

In this paper we propose a velocity-dependent 
Eddington factor /(t,/3) for a spherically symmetric case, 
and solve the fully special relativistic spherical outflows 
driven by radiation pressure using such /(r, /?). 

In the next section we propose a variable Eddington 
factor for spherical relativistic radiative flows. In section 
3 we describe the basic equations for relativistic outflow 
driven by radiation pressure under the spherical symme- 
try. In section 4 we show our numerical results of the 
radiative ffow. The final section is devoted to concluding 
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remarks. 

2. Variable Eddington Facor 

In this section we propose and explain an optical depth 
and velocity-dependent variable Eddington factor, which 
enables us to treat a problem of relativistic radiation hy- 
drodynamics in a spherically symmetric case. 

2.1. Traditional Eddington Factor 

We solve the radiation hydrodynamic problem scmi- 
analytically using moment equations. Then, the 
Eddington approximation is generally used to close mo- 
ment equations. The usual Eddington approximation is 
defined in the comoving frame as 

^0 = Ieo (1) 

where Pq is the radiation pressure and Eq is the energy 
density both measured in the comoving frame. This usual 
Eddington approximation is axiomatic when the radiation 
field is isotropic. Such a situation can be satisfied, when 
the atmosphere is sufficiently optically thick, or when the 
gas is optically thin with the uniform radiation field. As is 
well-known, however, this usual Eddington approximation 
does not hold when the radiation field becomes anisotropic 
in such a case that there is a transition from optically thick 
to thin states. 

When we examine the gas accelerated up to the rela- 
tivistic speed by radiation pressure, as a clue of mecha- 
nism for jets in active galactic nuclei and microquasars, it 
is necessary to consider a sudden change of optical depth 
and the steep velocity gradient. In order to obtain the 
terminal speed of the radiatively-driven relativistic out- 
flow, we have to investigate the flow down to the optically 
thin state. In addition, the usual Eddington approxima- 
tion would be violated, when the gas is accelerated to the 
relativistic speed with steep velocity gradient. This is be- 
cause the radiation field may become anisotropic, even 
in the comoving frame, due to relativistic aberration and 
redshift. 

2.2. Optical- Depth Dependent Factor 

In this subsection, we discuss about a better way of 
dealing with radiation field which is anisotropic. When 
there is a transition from optically thick to thin regimes, 
for a spherically symmetric case Tamazawa et al. (1975) 
set the Eddington approximation as 



understood as follows. The photon mean-free path £ is on 
the order of 



Po = .fEo 



(2) 



where / is a variable Eddington factor, and they proposed 
the relation that satisfied the physical condition from op- 
tically thick to thin regimes by 

/M^5^, (3) 

where T is the optical depth. 

This factor becomes 1/3 in an optically thick region 
while becomes unity in an optically thin region. This is 



e^i/{Kp) 



(4) 



where k is the opacity and p is the gas density. When the 
gas density is large and the medium is sufRciently thick, 
the mean free path becomes small and the radiation field 
is locally seen to be isotropic. While, around the surface of 
the atmosphere or in a spherically expanding flow, the gas 
density becomes small and the mean free path lengthens 
more and more toward the direction of the density gra- 
dient; then the radiation field becomes locally seen to be 
anisotropic. In such a transition region, the relation be- 
tween the radiation pressure and radiation energy would 
change in each direction. When the optical depth becomes 
0, for an outward direction, the radiation pressure is equal 
to the radiation energy. However, in the case of relativistic 
outflow, analytic method can not be calculated until the 
speed of light with even this factor due to the singularity. 

2.3. Velocity Dependent Factor 

Next, wc consider the case where the gas interacting 
with photon is accelerated to the relativistic speed. When 
there is a large velocity gradient, the photon mean- free 
path becomes longer than that without the velocity gra- 
dient. In such a case, the usual Eddington approximation 
would be violated again. For instant, in the relativistic 
flow with a velocity gradient dv/dr, where v is the flow 
velocity and r the radius, the velocity increase at a dis- 
tance of the mean free path i becomes 

dv 1 dv dv ,^, 

Av = e— = ——^—. (5) 

dr Kp dr dr 

In order for the radiation fields to be isotropic in the co- 
moving frame, this velocity increase should be sufficiently 
smaller than the speed of light; dv/dr ^ v/t <^ c. If the 
velocity difference becomes very large when the velocity 
itself is very high and/or the optical depth is small, the 
usual Eddington approximation in the comoving frame 
would be violated. Such a situation can occur for a rel- 
ativistic outflow. If the velocity difference is large at a 
distance of the mean free path, the relativistic effect, such 
as a Doppler effect and aberration, becomes important, 
and the radiation field is seen to be anisotropic. 

For a relativistic flow with a velocity gradient, a 
velocity-dependent variable Eddington factor was pro- 
posed (Fukue 2006): 



m-l^lp. 



(6) 



where f3 = v/c. This factor is applied to the plane-parallel 
case. As a good news by the usage of this velocity- 
dependent factor, we can avoid critical points that always 
appear in the moment equations under special relativity. 
In this point, it is indicated that the velocity-dependent 
variable Eddington factor in the relativistic flow could be 
reasonable mathematically as well as physically. However, 
in the spherical case we also have to consider about the 
effect of the optical depth against to the plane-parallel 
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case which does not include the effect of the optical depth 
through velocity of gas. 

2.^. Optical-Depth and Velocity Dependent Factor 

Now, we consider the case of a relativistic spherical flow. 
The Eddington factor depends on the optical depth for a 
spherical atmosphere, while it depends on the flow ve- 
locity for a relativistic flow. In the spherically symmet- 
ric relativistic flow, there exist a dilution effect due to a 
spherical expansion and that due to a relativistic expan- 
sion. Hence, we suppose that the Eddington factor could 
depend on both the optical depth and the flow velocity. 
The minimum requirements for such a variable Eddington 
factor are (i) it approaches 1/3 in a sufficiently thick, low 
velocity regime, (ii) it becomes unity in an optically thin 
regime, and (iii) it does also become unity in the relativis- 
tic regime at a speed on the order of the speed of light. 
Additional conditions are (iv) it reduces to the factor of 
Tamazawa et al. (1975) in a static limit, and (v) it is 
simple. 

Although there may be many possible factors, in the 
present paper we propose the following one. 



fir,P) 



7(l + /?) + T 
7(l-H/3)-f-3r' 



(7) 



where r is the optical depth, /3 is the normalized flow speed 
{j3 — v/c), and 7 is the Lorentz factor [7= l/-\/l — (w/c)^]. 
This form was born as follows. It is shown that the mean 
free path i of photons in the inertial frame lengthens 
than that io in the comoving frame by a relativistic ef- 
fect (Abramowicz et al. 1991) a.s i = £o/[7(l — /3cos9)] = 
^o7(l + /3)- By considering this, we replaced the optical 
depth of Tamazawa et al (1975) by T/[7(l-f /?)] for the 
outward moving flow. 

Figure 1 shows the behavior of the present variable 
Eddington factor (7). A dashed curve is the varial- 
blc Eddington factor by Tamazawa et al. (1975), while 
other curves are the present case for several values of 
the flow speed. As is seen in figure 1, the present vari- 
able Eddington factor becomes unity as the flow speed 
approaches the speed of light. 

Using these variable Eddington factors, we can calculate 
the spherically symmetric relativistic flow, continuously 
from low speed to relativistic regimes. In the next section, 
we solve the relativistic moment equations with a variable 
Eddington factor for the relativistic spherically symmetric 
case. 

3. Basic Equations 

In this paper, it is treated a simple one-dimensional 
radiation flow in what follows; i.e., we consider the spher- 
ical case in the radial direction. The radiative energy is 
transported in the radial direction, and the gas itself also 
moves in the radial direction by the action of radiation 
pressure. For simplicity, the radiation field is sufficiently 
intense that both the gravitational field, e.g., of the cen- 
tral object, and the gas pressure and the internal heating 
are ignored in this paper. As for the order of the flow vc- 
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Fig. 1. Optical deptli and velocity dependent variable 
Eddington factor. The dashed curve shows an optical depth 
dependent variable Eddington factor, /(r) = (1 + t)/(1-|-3t). 
The dotted lines denote the velocity dependent variable 
Eddington factor, /(/3) = 1/3 -I- (2/3)/3, where the veloc- 
ity is 13 = 0.3, 0.5, 0.7, 0.9, 0.99, 0.999, 0.99999 from bot- 
tom to top. The solid curves represent an optical 
depth and velocity dependent variable Eddington factor, 
/(r,/3) = {7(1-1-/3) -I- t}/{7(1 + /3) + 3t}, where the velocity 
is /3 = 0.3,0.5,0.7,0.9,0.99,0.999,0.99999 from bottom-left to 
top-right. 

locity V, we consider the fully special relativistic regime, 
where the all terms are retained. Importance to retain 
the higher order of velocity is described in Yin and Miller 
(1995). Under these assumptions, the radiation hydrody- 
namic equations for steady radial (r) flows arc described 
as follows (Kato et al. 1998; cf. Fukue 2006 for a plane- 
parallel case). 

The continuity equation is 



Anr^pcu — M {— const.). 



(8) 



where p is the proper gas density, u the radial four velocity, 
M the mass-loss rate, and c the speed of light. The four 
velocity u is related to the proper three velocity v hy u = 
7w/c. 

The equation of motion is 

c^u'^ = l^^lUls^ [F7(l + 2u') - c{E + P)^^u] , (9) 
dr c 

where Kabs and Ksca are the absorption and scattering 
opacities (gray), which relate to bremsstrahlung, photo- 
nionization and electron scattering. We define them in 
the comoving frame. Moreover, E is the radiation energy 
density, F the radiative flux, and P the radiation pressure 
observed in the inertial frame. In the no-gas pressure ap- 
proximation and without heating, the energy equation is 
reduced to a radiative equilibrium relation. 



= j - CKabs£'7 - CKabs^lt + 2Kabs-F7", 



(10) 



where j is the emissivity defined in the comoving frame. 
In this equation (10), the third and fourth terms on the 
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right-hand side appear in the relativistic regime. 

For radiation fields, the zcroth-nioment equation be- 
comes 

+K^bsFu/j~K^c^F{l + v^/c^)ju] . (11) 
The first-moment equation is 
dP 1 P7 

-J- = (3P - £■) H [ju/"f - Kabs^ + CKiihsPu/j 

ar r c 

-Ksca-F(l + 2U^) + CK^UE + P)iu] . (12) 

Although the first term on the right-hand side of equa- 
tion (12) disappears for the closure relation such as an 
Eddington approximation in the optically thick limit, this 
term remains in the equation because of the modification 
of the Eddington approximation in this paper. It should 
be noted that this equation (12) is reduced to that by 
Ruggles and Bath (1979) in the lower approximation of 

{v/cy. ^ 

Here, in order to close moment equations for radiation 
fields, we adopt a velocity-dependent or optical depth and 
velocity-dependent variable Eddington approximation (2). 
If we adopt this form (2) as the closure relation in the 
comoving frame, the transformed closure relation in the 
inertial frame is 

cP (1 + u^ - fu^) - cE {fj^ - u^) + 2Fju (1 - /) , (13) 

or equivalently, 

cP (1 - //32) ^cE{f- P^) + 2F/3 (1 - /) . (14) 

Above closure relation gives the relation among radiation 
pressure, energy and flux. Relations among E, F, and P, 
which depend on velocity, are important relations, and it 
is a point on using the modified closure relation. 

Eliminating j with the help of equations (10) and using 
continuity equation (8), equations (9), (11) and (12) are 
rearranged as 



dr 






X [F{l + 2u^)~c{E + P)-fu]. 



dr 



{r'^F) = r'^pu{K^bs + Ksca) 



X [c{E + P)-fu - F{1 + 2u^)] , 

(3P - E) + (Kabs + Ksca) 

r 

xp^ [c{E + P)ju - F(l + 2u^)] 



(15) 



(16) 



(17) 



The integration of the sum of equations (15) and (16) 
yields the energy flux conservation along the flow, 

c^Mj + L = c^M + La (= const.), (18) 

where L (= Airr'^F) is the luminosity. The initial condi- 
tions are given as u = 0, L = Lq, P — Pq, and r — tq at 
T = To . The subscript zero denotes the values at the flow 
base of T = tq. On the basis of above the basic equations 
are the equation of motion (15), the mass flux (8), the 



momentum flux (17), the energy flux (18) and the closure 
relation (14) at this stage. 

Here, we define new variables for convinient calcula- 
tions: Q = Attt^cP for radiation pressure and D = Aur^cE 
for radiation energy. Substituting these variables into 
equations (15), (17), and (18), with the help of equation 
(14), we obtain 



dr 



dQ 
dr 



(Kabs + Ksca)P7 

{f + /3^)L^/3Q{l + f) 

-('«abs + Ksca)/37 



(19) 



(/ + /?2)L-^Q(1 + /) 



l (l-/) (l + /3^)Q-2/3£(l-/) 
r 



(20) 
(21) 



In addition, we regard the optical depth t as 

dT=- (/tabs + Ksca)pdr, (22) 

and the mass flux (8), the momentum (19), the first mo- 
ment (20), and the energy flux (21) are rewritten as 

dr 47rr^C7/3 



dr 

dr 
dQ 
dr 



(/32+/)L-(l + /)/3Q 



-7- 



(23) 



(24) 



47rrc7/3 {1 -/)[{! + f3^)Q ~2f3L] 

(^abs + Ksca)M f - P'^ 

{f + P')L~PQil + f) 



Mc^-f + L = ]\Ic^ + Lo. (26) 

In order to transform them into dimentionless forms, 
the radius r is normarized by the Schwarzschild radius 
7'g (~ 2GAf/c^), the mass-loss rate M is normalized by 
-Le/c^, and the pressure Q and luminosity L are normal- 
ized by the Eddington luminosity Le [= 47rcG'M/(Kabs + 
'^sca)]- They can be rewrirren as 

df _ 2f27/3 
d^^ " 



Mj 



,dp 
dr 
dQ 
dr 



-1 



M 

{p^ + f)L-{l + f)PQ 



(27) 
(28) 

(29) 

(30) 

At this "first" normalization stage, we briefly comment 
the boundary conditions on the present case. Moment 
equations are to be solved as a two-point boundary value 
problem, as is well known. That is, at the base, flow deep 
inside the atmosphere, several conditions are imposed on 



2nP{l-f)[{l + p^)Q~2pL] 
{f + P^)L~PQ{l + f) 



-^1- 



f-P^ 



Mj + L = M + Lq. 
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the physical quantities for radiation fields, whereas, at the 
surface of the atmosphere, some relation generally holds 
on the radiative moments with or without the external 
irradiation. In the present radiative flow, we give the 
boundary conditions tq, /3(= 0), Qq (or Pq), and Lq at 
the flow base of the optical depth tq. In addition, there 
exists some relation for the moment Q and L (Fukue 2006) 
at the flow top of the optical depth of t = 0. Then, the 

mass-loss rate M should be determined as an eigcn value 
by the boundary condition at the flow top. Although it 
indeed be possible we do the "second" normalization be- 
low. 

In the present treatment, we only consider the radiation 
field without gravitational field add up to nothing char- 
acterisitic scale expected for mass-loss rate; i.e. mass-loss 
rate itself can be absorbed in the normalization unit. We 
further renormalize the variables by f = r/M, L = L/M, 

Q = Q/M to yield 

(IT 

= -2f^/?, (31) 



7 



dr 

dr 
dQ 



'1- 



(/32 + /)L-(l + /)/?Q 



d? = "^''^^ 



-/)[(! 4 



/32)Q_2/3L] 



(./ + /3')i- 



-/3Q(1 + ./) 



/-/32 



7 + i = 1 + io 



(32) 



(33) 
(34) 



At this "second" normalization stage, the mass-loss rate 
apparently disappears in the basic equations and it seems 
unnecessary the boundary condition at the flow top to 
determine the mass-loss rate. Thus, we solve equations 
(31)-(34) for a suitable form of variable Eddington factors 
/(t,/3). 

4. Results and Discussion 

In this section we briefly show a typical example for 
the rclativistic spherical flow using the present variable 
Eddington factor, and discuss and compare several forms 
of variable Eddinton factors. 

J^.l. Typical Example with Fastest Terminal Velocity 

Wc first show a typical example of the rclativistic spher- 
ical flow, after solving the special rclativistic radiation hy- 
drodynamic equations, using the present proposed factor 

(7). 

Among various combinations of parameters, we find the 
case of fastest terminal velocity for the initial condition at 
the flow base: Lq = 1, Qq = 0.99, fo = 1, /3o = 0, and tq = f . 
In this case, the terminal speed becomes 0.68c. The re- 
sult is shown in figure 2. As is seen in figure 2, the gas is 
accelerated as the luminosity decreases; i.e., the radiation 
energy is converted to the bulk motion in such a rclativis- 
tic regime. The gas is accelerated at around the flow top 
of f ~ blrgC^ / Le], where the optical depth vanishes. It is 
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Fig. 2. Flow three velocity v (thick solid curve), radiative 
luminosity L (thin solid curve), radiation pressure Q (chain- 
dotted curve), and radial distance r (doted curve) as a func- 
tion of the optical depth r for ro = 1 at the flow base of tq = 1. 
Other parameters are Lq = 1 and Qq = 0.99. 

stressed that there does not appear pathological critical 
points inherent in the usual Eddington factor of 1/3. 

Although we can flnd the rclativistic flow beyond the 
critical points, it is difficult to obtain the solutions with 
terminal velocities of '^ c. There are several reasons. 

The first is the restriction from energy conservation 
(34). Because the total energy is constant, even if all of 
the radiation energy is converted to the bulk energy, the 
terminal Lorentz factor 700 is restricted as 



7oo < 1+io- 



(35) 



In the case of figure 2, the terminal speed is smaller than 
this absolute limit. 

The second reason is the existence of radiation drag. On 
the right-hand side of equation (32), the term of L is the 
radiative acceleration, while the term related to Q means 
the radiation drag force, which is approximately propor- 
tional to the fiow speed. The terminal speed generally 
becomes high for large luminosities. At the same time, 
however, the radiation drag force becomes important as 
the flow speed is high. In the optically thin regime, radi- 
ation drag becomes important. 

The third is the dilution (curvature) effect of the spher- 
ical flow, which does not exist in the plane-parallel case 
(Fukue 2006). For a simple discussion, we assume that 
the fiow speed is constant with the terminal value of /3oo 
and 7oo. In such a case, the continuity equation (31) is 
integrated from the flow base of tq to top of Vqc as 

11 

= 27oo/3ooTo. (36) 

ro roo 

Or, there is a restriction of 



loo Poo < 



2roTo 
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for finite rg . The terminal speed in figure 2 is on the order 
of this dilution (curvature) limit. 

4-. 2. Comparison with Various Factors 

In this subsection we compare the results for various 
variable Eddington factors. 

In figures 3-5 the velocity, the luminosity, and the ra- 
dius are shown respectively as a function of the optical 
depth for parameters of Lq = 1, Qq = 1, and tq = 1. The 
thin solid curve is results for usual Eddington factor which 
is constant. The chain-dotted curve is the one for the 
factor of Tamazawa et al.(1975) which depends on opti- 
cal depth. The dashed curve is the one for the factor is 
Fukue(2006) which depends on velocity. The thick solid 
curve is the one for the factor of proposed factor in this 
paper which depends on velocity and optical depth. 

As seen in figure 3, the result for usual Eddington factor 
shows an acceleration of the flow is not enough due to the 
singular or radiation drag. On the other hand, in the case 
of the other factors it can be also shown that the terminal 
velocity of gas becomes large, even for the same param- 
eters. In this case there is a difference of around 15% in 
each factor. Although it seems to close these results for 
Tamazawa's factor and present proposed one, the latter is 
physically acceptable as discussed in section 2. 

As seen in figure 4, the luminosity is converted to the 
bulk motion. That is, the luminosity decrease for the 
usual factor is small, whereas that for other factors is large 
up to 10% ~ 15%. Such a luminosity change may be a 
clue to discriminate the various Eddington factors. 

Finally, figure 5 shows the radius change of the expand- 
ing photosphere. The terminal speed is large as the accel- 
erating distance becomes large. In the present example, 
the radius is at most f ~ ^rgC^ /Le]- For large optical 
depth, this radius would be large, and the terminal speed 
would also become large. These results show that gas ac- 
celerate as a stretch in the vicinity of the center of compact 
objects. 

5. Concluding Remarks 

In the present paper, we examine the relativistic ra- 
diation fiow in the spherically symmetric case with the 
velocity- and optical depth- dependent variable Eddington 
factors within the framework of special relativity. We 
showed that in the relativisic spherical fiow the Eddington 
factor is no longer constant, but depends on the velocity 
as well as the optical depth. In particular, when the gas 
is accelerated up to the relativistic speed, there exists a 
strong velocity gradient, and the velocity dependence of 
the Eddington factor becomes important. In addition, 
such a variable factor can avoid the pathological singu- 
larity in the moment equations. We emphasize that we 
should use such a generalized Eddington factor to treat 
the relativistic radiation hydrodynamics under the mo- 
ment formalism. 

We can find several solutions for the relativistic spheri- 
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Fig. 3. Velocity as a function of tlio optical depth for pa- 
rameters of Lo = 1, Qo = 1, and rg = 1. The thin solid 
curve is for the case of / = 1/3, the chain-dotted one for 

/(t) = (1 + t)/(1 + 3t), the dashed one for /(/3) = (1 + 2/3)/3, 
the thick solid one for /(t,/3) = {7(1-1- /3) + t}/{7(1 + /3)-|-3t}. 
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Fig. 4. Luminosity as a function of the optical depth for 
parameters of Lq = 1, Qo = 1, and ro = 1. The thin solid 
curve is for the case of / = 1/3, the chain-dotted one for 
/{T) = {l + T)/(l-|-3r), the dashed one for /{/?) = (1 + 2/3)/3, 
the thick solid one for /(t,/3) = {7(1-1-/3) + r}/{7(l-f/3)-|-3T}. 
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Spherical Rclativistic Radiation Flows 




Fig. 5. Radius as a function of the optical depth for pa- 
rameters of Lo = 1) Qo = li and ro = 1. The thin sohd 
curve is for the case of / = 1/3, the chain-dotted one for 
/(r) = (l-|-r)/(l-|-3T), the dashed one for /(/3) = (l-|-2/3)/3, 
the thick solid one for /(t,/3) = {7(l-|-/3)-|-r}/{7(l-|-/3)-|-3T}. 

cal flow. The results, however, are slightly different for the 
Eddington factor adopted. In order to determine the pre- 
cise form of the Eddington factor, we must solve the rel- 
ativistic transfer equation rigorously. However, the func- 
tional form of the variable Eddington factor is usuful for 
the study of the relativistic jets, black-hole winds, and the 
gamma-ray bursts. 

It should be commented on the current works on the 
related topics. Current works are divided mainly into 
two categories, as refered in the introduction. One type 
solved the relativistic radiation hydrodynamical equations 
under the diffusion approximation (e.g., Ruggles, Bath 
1979; Quinn, Paczyiiski 1985; Paczyhski, Proszyhski 1986; 
TuroUa et al. 1986; Paczyhski 1990; Nobih et al. 1994). 
In these current works the flow is restricted in the sub- 
relativistic region on the order of ~ 0.1 c. However, the 
diffusion approximatioon may be valid only in the suf- 
ficiently optically thick regime, and further, there is no 
justification that the diffusion approximation can be used 
in the relativistic regime, since there exists a causality 
problem. Another type examined the pathological behav- 
ior of the traditional moment formalism in the relativis- 
tic regime (e.g., TuroUa, Nobih 1988; Nobih et al. 1991; 
TuroUa et al. 1995; Dullemond 1999; Fukue 2005), which 
is one of the motivation of the present study. However, 
there is no proposal to use a variable Eddington factor, 
which depends on the fiow velocity as well as the optical 
depth, in order to solve the moment equations in the rela- 
tivistic regime in the spherically symmetric case. We thus 
tried to solve the relativistic moment equations with an 
approximate form of the variable Eddington factor. 

In this paper, we considered only the one-dimensinal 
case without gravity under special relativity. In order to 
clarify the physics of rclativistic jets around a black hole. 



we must treat the problem within the framework of gen- 
eral relativity. Such a case is a next work. 

This work has been supported in part by a Grant-in- 
Aid for Scientific Research (18540240 JF) of the Ministry 
of Education, Culture, Sports, Science and Technology. 
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